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1. Introduction 



Quantum Chromodynamics is a remarkable theory. It is a convincing practical example of the 
triumph of the quantum field theory. Asymptotic freedom allows QCD to be consistent down to 
arbitrary short distance scale, enabling us to define the theory completely in terms of the funda- 
mental microscopic degrees of freedom - quarks and gluons. This fundamental definition is very 
simple, yet the theory describes a wide range of phenomena - from the mass spectrum of hadrons 
to deep-inelastic processes. As such, QCD should also possess well defined thermodynamic prop- 
erties. The knowledge of QCD thermodynamics is essential for the understanding of such natural 
phenomena as compact stars and laboratory experiments involving relativistic heavy-ion collisions. 

Full analytical treatment of QCD is very difficult because, neglecting quark masses, this theory 
has no numerically small fundamental parameters. The only independent intrinsic scale in this the- 
ory is the dynamically generated confinement scale Aqcd ~ lfm -1 . In certain limits, in particular, 
for large values of the external thermodynamic parameters temperature T and/or baryo-chemical 
potential /ig, when thermodynamics is dominated by short-distance QCD dynamics, the theory can 
be studied analytically, due to the asymptotic freedom. But, as seen below, the most interesting 
experimental region of parameters T and jU# is that of order Aqcd- 

The above makes first principle lattice approaches, which do not rely on small parameter ex- 
pansions, an invaluable and the most powerful tool in studying QCD thermodynamics. In addition, 
the domain where all relevant time/distance (or energy /momentum) scales are similar is especially 
suited for a lattice study (accommodating a wide scale window would require a correspondingly 
large lattice). 

The full potential of lattice methods is close to being realized as far as the study of QCD at 
IXb = is concerned. The main practical problems in this regime - accommodating sufficiently light 
quarks and approaching the continuum limit - are being methodically and successfully addressed 
through the use of improved lattice discretization schemes, as well as advances in algorithm and 
hardware technology. 

The status of thermodynamics of QCD at non-zero fls is different. The main impediment to 
lattice simulations is the notorious sign problem, discussed in Section 3.1. No method devised so 
far is known, or expected, to converge to the correct physical result as the infinite volume limit is ap- 
proached at fixed Hb 7^ 0. However, since the most interesting structure of the QCD phase diagram 
(phase transitions and critical points) lie at nonzero any progress in this direction is especially 
valuable. Existing lattice methods generically rely on clever extrapolations from = 0. These 
techniques yield interesting results in the regime of small, but already experimentally relevant jx B - 

A contemporary view of the QCD phase diagram is shown in Fig. ||[ It is a compilation of a 
body of results from model calculations, empirical nuclear physics, as well as first principle lattice 
QCD calculations and perturbative calculations in asymptotic regimes. 

Several reviews in these proceedings, in addition to original contributions, are devoted to re- 
cent progress in lattice studies of QCD thermodynamics. Ref. [|l|] reviews thermodynamics of 
QCD at Hb = 0. Ref. [^] discusses lattice results at small Hb- In addition, Ref.[Q] describes resent 
progress in uncovering phase structure of QCD at large relevant to the physics of compact stars, 
and outlines targets of opportunity for potential future lattice studies in this domain. 
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This report provides an overview of the structure of the QCD phase diagram based on available 
theoretical (lattice and model calculations) and phenomenological input. Some of the recent lattice 
results reported separately in this volume are also briefly discussed in Section Q 

2. The phase diagram 

Thermodynamic properties of a system are most readily expressed in terms of a phase diagram 
in the space of thermodynamic parameters - in the case of QCD - as a T/ig phase diagram. Each 
point on the diagram corresponds to a stable thermodynamic state, characterized by various ther- 
modynamic functions, such as, e.g., pressure, baryon density, etc (as well as kinetic coefficients, 
e.g., diffusion or viscosity coefficients, or other properties of various correlation functions). 

Static thermodynamic quantities can be derived from the partition function - a Gibbs sum 
over eigenstates of QCD Hamiltonian, which can be alternatively expressed as a path integral in 
Euclidean space: 

Z{T^ B ) = L^p{~ Ea ~t BB " } = / ^0W)exp{-S*} (2-D 

where a labels states with energy E a and baryon number B a . The path integral is over color gauge 
(gluon) fields periodic in Euclidean time with period l/T, and quarks fields, antiperiodic with 
the same period. The Euclidean action given by 

S E = Sym + £ f d 4 x q^Oq (2.2) 

q=u,d,s 

where Sym[^] * s trie SU(3) Yang-Mills action and, in the chiral Weyl basis, the Dirac spinors and 
matrix are 

1=( qL ) and B=(°- D 7* (2.3) 
where cfy = (1, ia), = — iA^, and \i q = /^b/3. 

2.1 Chiral symmetry argument 

In the chiral limit - the idealized limit when 2 lightest quarks, u and d, are taken to be massless, 
the Lagrangian of QCD acquires chiral symmetry SU(2)ixSU(2)s, corresponding to SU(2) flavor 
rotations of (mz,,^l) and (mr,^) doublets independently. The ground state of QCD breaks the chiral 
symmetry spontaneously locking SU(2)^ and SU(2)« rotations into a single vector-like SU(2)y 
(isospin) symmetry and generating 3 massless Goldstone pseudoscalar bosons - the pions. The 
breaking of the chiral symmetry is a non-perturbative phenomenon. 

At sufficiently high temperature T ^> Aqcd> due to the asymptotic freedom of QCD, pertur- 
bation theory around the approximation of the gas of free quarks and gluons (quark-gluon plasma 
- QGP) should become applicable. In this regime chiral symmetry is not broken. Thus we must 
expect a transition from a broken chiral symmetry vacuum state to a chirally symmetric equilibrium 
state at some temperature T c ~ Aqcd- The transition is akin to the Curie point in a ferromagnet - 
where the rotational 0(3) symmetry is restored by thermal fluctuations (chiral 0(4)=SU(2)xSU(2) 
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Figure 1: Phase diagram of QCD with massless quarks dictated by the chiral symmetry argument. The order 
of the transition (solid magenta line) is not determined by this simplest argument. 

symmetry in QCD). Thermodynamic functions of QCD must be singular at the transition point - 
as always when the transition separates thermodynamic states of different global symmetry. 

Thus, the region of broken chiral symmetry on the T/ig phase diagram must be separated from 
the region of the restored symmetry by a closed boundary as shown in Fig. |I[ 

2.2 Pisarski-Wilczek argument 

The chiral symmetry argument alone is not sufficient to determine the order of the temperature 
driven chiral symmetry restoration transition. A more elaborate argument, based on universality, 
advanced by Pisarski and Wilczek [|J] asserts that the transition cannot be of second order for three 
massless quarks. 

In a simplified form, the logic of Ref.[Q] is as follows. Let us assume that the transition is of 
the second order. Then the critical behavior of the system (long-distance behavior of correlation 
functions, singular contributions to thermodynamic functions, etc.) is determined by the long- 
wavelength modes which, in the case of the second order transition in a theory with Nf light quarks, 
are the Nf — 1 pions of the spontaneously broken S\J(Nf)A axial flavor symmetry plus the critical 
mode - the magnitude of the chiral condensate <T ~ qq. 

Universality implies that the critical behavior is the same as in any local theory in 3 dimensions 
with the same global symmetry breaking pattern and the same set of critical modes. In our case, 
a representative example of the universality class is an SU(7Vf)x SU(Aff) sigma model of an Nf x 
Nf matrix- valued field <1>. It turns out, that for Nf = 3, the model cannot be critical: there is a 
relevant operator cubic in the order parameter field, det<l>, which always destabilizes the symmetric 
minimum of the effective potential for <1> via a first order transition, before the curvature of the 
minimum vanishes (i.e., before criticality is reached). Hence, QCD with Nf = 3 massless quark 
flavors must undergo a. first order chiral restoration transition. 

2.3 Nf = 2 chiral limit and tricritical point 

For two massless quarks the transition can be either second or first order. As lattice and model 
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2 massless quarks 




1 \i B , GeV 



Figure 2: The order of the chiral symmetry restoration transition for 2 massless and one massive quarks. 
The dashed line (red) is the second order transition, the solid line (blue) is the first order transition. In the 
low T region: chiral symmetry is broken in nuclear matter. Details of the phase structure at high /Ib are 
omitted. 



calculations show, both possibilities are realized depending on the value of the strange quark mass 
m s and/or the baryo-chemical potential Hb- 

The point on the chiral phase transition line where the transition changes order is called tri- 
critical point, see Fig. ||. The location of this point is one of the unknowns of the QCD phase 
diagram with 2 massless quarks. In fact, even the order of the transition at }i B = 0, which many 
older and recent studies suggest is of the second order (as shown in Fig. |2|) is still being questioned 
(see review by Heller in this volume [{[]]). 

Neither can it be claimed reliably (model or assumption independently) that the transition, if it 
begins as a 2nd order at jj, B = 0, changes to first order. However, numerous model calculations show 



this is the case (Section 3.4). Lattice calculations also support such a picture. Recent advances 
in the understanding of QCD at low T and large Hb, reviewed in [Qj, also point at a first order 
transition (at low-7 1 , high-/^) from nuclear matter to color-superconducting quark matter phase. 
Fig. |2| reflects this consensus. 

At low temperature, nuclear matter (which is expected to be still bound in the chiral limit) 
should be placed on the broken symmetry side of the chiral transition line as shown in Fig. fj. 

2.4 Physical quark masses and crossover 

When the up and down quark masses are set to their observed finite values, the diagram as- 
sumes the shape sketched in Fig. [| The second order transition line (where there was one) is 
replaced by a crossover - the criticality needed for the second order transition in Fig. || requires 
tuning chiral symmetry breaking parameters (quark masses) to zero. In the absence of the exact 
chiral symmetry (broken by quark masses) the transition from low- to high-temperature phases of 
QCD need not proceed through a singularity. Lattice simulations do indeed show that the transition 
is a crossover for jJ. B =0 (most recently and decisively Ref.[Bt|, see also Ref.[]T]] for a review). 1 



'This fact is technically easier to establish than the order of the transition in the chiral limit - taking the chiral limit 
is an added difficulty. 
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Figure 3: The contemporary view of the QCD phase diagram - a semiquantitative sketch. 

This transitional crossover region is notoriously difficult to describe or model analytically 
- description in terms of the hadronic degrees of freedom (resonance gas) breaks down as one 
approaches crossover temperature (often called T c ), and the dual description in terms of weakly 
interacting quarks and gluons does not become valid until much higher temperatures. Recent ter- 
minology for the QCD state near the crossover (T ~ (1 — 2)T C ) is strongly coupled quark-gluon 
plasma (sQGP). 

Transport properties of sQGP have attracted considerable attention. For example, generally, 
the shear viscosity r\ is a decreasing function of the coupling strength. The dimensionless ratio of 
r\/fi to the entropy density s tends to infinity asymptotically far on either side of the crossover - in 
dilute hadron gas (T — > 0) and in asymptotically free QGP (T — > oo). Near the crossover r\/s should 
thus be expected to reach a minimum [g]. The viscosity can be indirectly determined in heavy 
ion collisions by comparing hydrodynamic calculations to experimental data. Such comparison [|7|] 
indeed indicates that the viscosity (per entropy density) of this "crossover liquid" is relatively small, 
and plausibly is saturating the lower bound conjectured in 

2.5 Physical quark masses and the critical point 

The first order transition line is now ending at a point known as the QCD critical point or end 
point. 2 The end point of a first order line is a critical point of the second order. This is by far 
the most common critical phenomenon in condensed matter physics. Most liquids possess such 
a singularity, including water. The line which we know as the water boiling transition ends at 
pressure p = 218 atm and T = 374°C. Along this line the two coexisting phases (water and vapor) 
become less and less distinct as one approaches the end point (the density of water decreases and 
of vapor increases), resulting in a single phase at this point and beyond. 

In QCD the two coexisting phases are hadron gas (lower T), and quark-gluon plasma (higher 
T). What distinguishes the two phases? As in the case of water and vapor, the distinction is only 

2 The QCD critical point is sometimes also referred to as chiral critical point which sets it apart from another known 
(nuclear) critical point, the end-point of the transition separating nuclear liquid and gas phases (see Fig. ^|). This point 
occurs at much lower temperatures ff(lOMeV) set by the scale of the nuclear binding energies. 
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quantitative, and more obviously so as we approach the critical point. Since the chiral symmeUy is 
explicitly broken by quark masses, the two phases cannot be distinguished by realizations (broken 
vs restored) of any global symmetry. 3 

It is worth pointing out that beside the critical point, the phase diagram of QCD in Fig. || has 
other similarities with the phase diagram of water. A number of ordered quark matter phases must 
exist in the low-T high-/ig region, which are akin to many (more than 10) confirmed phases of 
ice. For asymptotically large QCD with 3 quark flavors must be in color-flavor locked (CFL) 
state g |. 



3. Locating the critical point: theory 

The critical point is a well-defined singularity on the phase diagram, and it appears as an attrac- 
tive theoretical, as well as experimental, target to shoot at. Theoretically, finding the coordinates 
(T,jx B ) of the critical point is a straightforwardly defined task. We need to calculate the partition 



function of QCD given by Eq. ( |2.1[ ) and find the singularity corresponding to the end of the first 
order transition line. But it is easier said than done. 

Of course, calculating such an infinitely dimensional integral analytically is beyond present 
reach. Numerical lattice Monte Carlo simulations is an obvious tool to choose for this task. At 
zero pLjj Monte Carlo method allows us to determine the equation of state of QCD as a function of 
T (and show that the transition is a crossover). However, at finite [ig the Nature guards its secrets 
better. 



3.1 Importance sampling and the sign problem 

The notorious sign problem has been known to lattice Monte Carlo experts since the early 
days of this field. Calculating the partition function using Monte Carlo method hinges on the 
fact that the exponent of the Euclidean action Se is a positive-definite function of its variables 
(values of the fields on the lattice). This allows one to limit calculation to a relatively small set of 
field configurations randomly picked with probability proportional to the value of exp(— Se)- The 
number of such configurations needed to achieve reasonable accuracy is vastly smaller than the total 
number of possible configurations. The latter is exponentially large in the size V of the system, or, 
the number of the degrees of freedom: exp (const -V). The method, also known as importance 
sampling, utilizes the fact that the vast majority of these configurations contribute a tiny fraction 
because of the exponential suppression by exp(— Se)- Only configurations with sizable exp(— Se) 
are important. 

In QCD with jig ^ the Monte Carlo action Smc (playing the role of Se) is complex. With 
Smc complex, how does one pick important configurations? A number of ways to circumvent 
the problem have been tried. For example, using the modulus of exp (—.Smc) as a measure of 
importance, or the value of exp(— Smc) at zero ^Ug, when it is still positive. Unfortunately, none 
of the methods can be expected to converge to correct result with the increasing lattice volume, 

3 Deconfinement, although a useful concept to discuss the transition from hadron to quark-gluon plasma, strictly 
speaking, does not provide a distinction between the phases. With quarks, even in vacuum (T = 0) the confining 
potential cannot rise infinitely - a quark-antiquark pair inserted into the color flux tube breaks it. The energy required to 
separate two test color charges from each other is finite if there are dynamical quarks. 
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unless this limit is not accompanied by an exponential exp (const -V) increase of the number of 
configurations, rendering Monte Carlo technique useless. 

3.2 The overlap problem 

To demonstrate the problem, consider the most straightforward attempt to circumvent it - 
reweighting. 4 We cannot obtain a correctly weighted sample of important configurations directly 
at jUg ^ 0. But we still can at = 0. So, we take the {J-b = sample and offset incorrect probability 
by multiplying the contribution of each configuration by a factor exp(+5£ |^ B= o — Se)- This is exact 
in the limit when the sample contains all possible configurations. 

The problem is in the size of the sample needed for a Monte Carlo computation as V — > °°. 
The method uses the fact that at finite volume V, even at = 0, the configurations important 
for Hb y^O pop up, but with a very small probability. This probability is exponentially small as 
volume V — ► °°: exp (—const -V). When we calculate the partition function the reweighting factor 
is correcting for that, and is therefore exponentially large (for the complex Smc, both he magnitude 
and the complex phase are). Fluctuations, or statistical noise, in the exponentially tiny number of 
the rare important configurations completely washes out the significance of the result. 

In layman's terms, imagine that we want to study ice, but can only run experiments at normal 
room temperature and pressure. Using the reweighting method is analogous to trying to glimpse 
the information by waiting for rare configurations when all the water molecules accidentally gather 
in one corner of the lab, forming a chunk of ice. The amount of time that this experiment would 
require is exponentially large as V — > °°. 

3.3 Complex determinant 

Why is the Monte Carlo action in QCD complex and what can be done about it? To see, 



integrate over the quark fields in ( 2. 1 ) explicitly and obtain 



Z= f @Ae- SYM Y\detB= j S>A e ^ MC (3.1) 

J q J 



where (as in Eq. (2.3) and using the property = —D) 




(3.2) 



For H q = each quark determinant in Eq. 

detB = det[(a-D)(a-D) T + m^] > 0. (3.3) 

The positivity (and even reality) is lost if \i q ^ 0. This is the sign problem. 

However, the following still holds (detD(ju^))* = detB(— ju*). This opens two possibilities 



for the measure in the Euclidean path integral (3.1) to remain positive for fi q ^ 0: (a) if \i q is 



imaginary; or (b) if there are two degenerate quarks, e.g., m u = and \i u = —jXd, which is what 
happens with the chemical potential pLj of isospin ^3, or in phase-quenched QCD. Both alternatives 
are being exploited to glimpse into the regime }Xb / 0, yet unaccessible to direct Monte Carlo. In 



particular, the recent results from the simulations at finite jij are reported in Ref.[l 1 1. Simulations 



at imaginary \Xb are discussed further in Section ^2 



For QCD at finite [1b this method is known as the "Glasgow method" (reviewed in Ref. ||10||). 
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Figure 4: Comparison of predictions for the location of the QCD critical point on the phase diagram. Black 
points are model predictions: NJLa89, NJLb89 - [jlf], C094 - ^ [wj, INJL98 - JIJ], RM98 - [jTfj], 
LSM01, NJL01 - 10, HB02 - CJT02 - [Q, 3NJL05 - @, PNJL06 - Green points are lattice 
predictions: LR01, LR04 - [g, LTE03 - [pj, LTE04 - pi. The two dashed lines are parabolas with 



slopes corresponding to lattice predictions of the slope dT jd[i\ of the transition line at fig — p4 
The red circles are locations of the freezeout points for heavy ion collisions at corresponding center of mass 
energies per nucleon (indicated by labels in GeV) - Section g. 



3.4 Predictions from models 



In the absence of a controllable (i.e., systematically improvable and converging in the V — > °° 
limit) method to simulate QCD at nonzero /ig, one turns to model calculations. Many such calcula- 
tions have been done [12, 13, 14], 15, T^, 17, 18, [19|, 20, 21]. Figure |] summarizes the results. One 
can see that the predictions vary wildly. An interesting point to keep in mind is that each of these 
models is tuned to reproduce vacuum, T = Hb = 0, phenomenology. Nevertheless, extrapolation to 
nonzero /ig is not constrained significantly by this. In a loose sense, most lattice methods (see next 
Section) can be also viewed as extrapolations from Hb = 0, albeit with reliable input from finite T. 



4. Lattice results on the critical point 

This section is devoted to brief (and necessarily incomplete) descriptions of currently devel- 
oped lattice methods for reaching out into the Tjj, B plane. The comments below are selective and 
are meant to complement the original contributions in this volume. For a more comprehensive 
description of these methods, as well as other methods not discussed here, the reader may consult 
the most up-to-date review of Schmidt in these proceedings [Q] as well as an earlier review by 
Philipsen p6[], both of which also contain further references to original papers. 

4.1 Re weigh ting 

The first lattice prediction for the location of the critical point was reported by Fodor and 
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Figure 5: Expected direction of motion of the critical point (left) as the quark masses are decreased as shown 
in the Columbia plot in the inset. 



Katz in Ref. [22]. The assumption is that, although the problem becomes exponentially difficult 
as V — > °°, in practice, one can get a sensible approximation at finite V. In addition, simulations 



at finite T might suffer lesser overlap problem because of large thermal fluctuations Q27|]. One can 
hope that if the critical point is at a small value of the volume V may not need to be too large to 



achieve a reasonable accuracy. In particular, numerical estimates show [ ]28| ] that the maximal value 
of /ig which one can reach within the same accuracy shrinks only as a power of l/V . 



The results of Ref. [ |22| ] are the most definitive and well-known, but they also attract the 
strongest of criticisms. The method of the Ref. [ |22| ] is based on computing the position of the 
zero of the partition function in the complex temperature plane and observing when (for which 
this zero crosses (and with its complex conjugate - pinches) the real axis. This determines the T 



and Hb coordinates of the critical point. However, as Ejiri points out in Ref. [29], once the fluctu- 
ations of the phase, argdetB, of the Dirac determinant are large, they cause fake zeros to appear. 



It is therefore alarming that, as Splittorff argues [ftOH , both points found in Ref. Q22Q (different m q 
and V) happen to lie on the critical line of the phase quenched QCD (| detB| instead of detB) - 
which is the line where fluctuations of argdetB do become large. In a related observation, Goiter- 
man et al pi], [32|] argue that the procedure of taking the fourth root [33] of the staggered fermions 



causes problems in a finite \1b calculation such as in Ref. [22] 



4.2 Imaginary [i B and Nf = 3 

By the universality argument of Section ^2|, the finite temperature transition is 1st order for 
m u = m d = m s = 0- By continuity, it must remain 1st order in a finite domain of the m s m uc i plane 
(taking m u = m c i = m u d) surrounding the origin - the plot of this domain is known as Columbia 
plot [|34|, [I]]. For physical quark masses and Hb = the temperature driven transition is a crossover, 
which means that the physical point is outside the 1st order domain in the m s m uc i plot. Reducing 
quark masses should pull the point into the 1st order domain. 

What happens on the T/ig phase diagram as the point m s m U( j is pulled towards and into the 
first order domain? The most straightforward expectation is that the first order line begins earlier, 
at lower i-e., the critical point is pulled towards the /j.b = axis, as shown in Fig. |5[ until 
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it disappears off the phase diagram altogether, and the whole transition line is of the 1st order. 5 
Furthermore, lattice calculations at /ig = show that real QCD is very near the 1st order domain 
boundary. That suggests the critical point is not too far off the T axis in the T/ig plane. 

What happens to the critical point when {m s ,m u d) is in the 1st order domain? It is still a 
singularity of the partition function as a function of Hb, but it moves out into the complex jU# 
plane. More precisely, it moves to imaginary }Xb axis. This remarkable fact allows one to observe 
the (complex descendant of) the critical point in a direct Monte Carlo simulation - since there is 



no sign problem for imaginary \Xb (Section 3.3). This observation is at the core of the method 



developed by de Forcrand and Philipsen [25, 35]. 

The success of the method crucially depends on the analyticity of the coordinate /ij of the 
critical point as a function of quark mass, e.g., m s around the point where \l\ = 0. The validity of 
this can be argued as follows. In the (T,HB,m s ) space the criticality is achieved (correlation length 
goes to infinity) when 2 conditions are satisfied: t{T,}Xg,m s ) = h(T ',/iJ,m iS .) = 0, i.e., there are two 
relevant operators in the universality class of the critical point and their coefficients, t and h, must 
be tuned to zero. The coefficients of these operators are analytic functions of the parameters. 6 
Furthermore, the analyticity in /ij (not just in jx B , which otherwise could cause a branching point 
at jiij = 0) follows from the — ► —jis symmetry of the QCD partition function. Solving the two 
conditions for T and /ij one finds the position of the critical point (T (m s ) , Hg(m s )) in terms of 
functions analytic in m s . 

de Forcrand and Philipsen determine the function /i|(m v ), or rather its inverse m s (jU.J), for 
/ij < and then analytically continue to real /i^. This way one could estimate the position of the 
critical point in the TjXs plane. 

It is puzzling that the slope of the function m s (Hjj) measured in this way [|25|, ^5|] appears 
negligible in lattice units and has a wrong sign 7 after a translation to physical units is applied. This 



leads the authors of Refs. [|25j, |35|] to suggest an unusual scenario: a new critical point is emerging 
on the phase diagram as the (m s ,m U( i) point is taken into the 1st order domain on the Columbia 
plot, dragging a new line of 1st order transitions into the T/ig plane. An unusual feature of such a 
point worth pointing out is the positioning of the 1st order line on the high temperature side of the 



critical point. As emphasized in Ref.[35], these results should still be subject to large uncertainties 
due to discretization and/or finite volume errors and more refined simulations are needed before 
physical conclusions can be drawn. 

4.3 Taylor expansion 

Taylor expansion in jU# is another method to circumvent the sign problem. Derivatives of 
pressure (or other thermodynamic quantities) are calculated at /j,b = and assembled into a Taylor 



series expansion to obtain dependence of that quantity on /ig [23, 24, 36] 

Consistent with the existence of the critical point at finite jj, B , there is a noticeable rise in the 
baryon number susceptibility Xb - see the peak on Fig. |[ Such a peak should be expected since the 



5 One can see how this scenario is realized on a 3-flavor NJL model in Ref,|pO[]. 

6 The non-analyticity characteristic of the critical behavior, arises due to non-analytic dependence of the correlation 
length, £ , on the values of the relevant parameters / and h. 
7 I.e., opposite to the sign implied by Fig. ^| 
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baryon number susceptibility diverges at the critical point. On the other hand, the isospin suscepti- 
bility should not diverge at the critical point, because the critical mode, a ~ qq is an isoscalar and 
cannot be excited by the operator of isospin /3 (while it is excitable by the operator of the baryon 
number). 



The authors of Ref. [ p7| ] caution against attributing the peak to the critical point. Their reser- 
vation is due to the fact that the low-J side of the peak is well described by a hadronic resonance 
gas model. Nevertheless, the agreement with the resonance gas does not necessarily mean that the 
rise of susceptibility cannot be due to the critical point of QCD. On the contrary, viewing resonance 
gas conceptually as a complementary (dual) description of QCD one must conclude that resonance 
description must reproduce the same thermodynamic functions as fundamental QCD description - 



including the critical point. Although the simple resonance gas model used in Ref. ||37|] does not 
describe the critical point itself, it still might be describing the onset of the critical behavior, just 
before the model breaks down. 8 

Furthermore, the resonance gas model of Ref.[]37|] does not describe the higher T side of the 
peak. The model must break down as the peak is approached from below, and is certainly not valid 
above the peak, where a different description must be used. At the same time, both sides of the 
peak receive a natural interpretation in terms of the proximity of the critical point. 

4.4 Radius of convergence of the Taylor expansion 

At a fixed temperature, the convergence radius of the Taylor expansion in Hb is limited by the 
nearest singularity in the complex plane of }Xb- Assuming that at the temperature Tg, at which the 
critical point (Te,He) occurs on the phase diagram, this critical point is the nearest singularity to 
jig = 0, one could use Taylor expansion to determine \i E [23, 24, [36|, 2|], if T E is known. 



Assuming that the radius of convergence }Xr can be approximated using the first few terms 
of the Taylor expansion, one can plot }Xr as a function of T. The main remaining problem is to 



8 A calculation illustrating this point has been reported earlier in Ref. Jig]: Improving the resonance gas description 
by a certain bootstrap procedure one obtains an equation of state which does have a critical point, similar to the van der 
Waals equation of state. 
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Figure 7: Refpq]: The location of complex plane singularities (cuts in the thermodynamic limit)) in the 
random matrix model at two representative values of T: T c and Tg. At Te the branching points pinch the real 
axis - this is where the critical point appears on the phase diagram. The trajectories of the branching points 
are shown by a dashed (red) curve. On the right the radius of convergence of the Taylor series, set by the 
distance from jig = to the branching point, is plotted as a function of T (the latter is along the ordinate to 
facilitate the comparison with QCD phase diagram). The critical point is shown. 



determine the value of Te i.e., to identify at which value of T the complex singularity reaches the 
real axis in the jUg plane. This question has been addressed using universality arguments, as well 
as an example random matrix calculation in Ref[38]. The trajectory of the complex singularities is 
illustrated in Fig. [7|. 

Two conclusions can be made: (i) The minimum value of the radius of convergence is not 
achieved at T = Te, but rather at a temperature close to the temperature T c of the chiral transition 
at Hb = 0. 9 (ii) At Te the function Hr(T) has a high order singularity. 

It is unlikely that such a weak singularity alone can be used to identify the value of Te- This 
suggests that one should attempt to extract more information from the Taylor series, for example, 
using the complex phase of the /Xg-plane singularity at given T . The critical point could then be 
located by the condition that this singularity is on the real axis. Such analysis would require observ- 
ing sign oscillations of the Taylor coefficients, and will require the knowledge of the coefficients 
up to an order higher than available to date. 



5. Scanning QCD phase diagram in heavy ion collisions 

Even though the exact location of the critical point is not known to us yet, the available the- 
oretical estimates suggest that the point is within the region of the phase diagram probed by the 
heavy-ion collision experiments. This raises the possibility to discover this point in such experi- 
ments [J3|]. 

It is known empirically that with increasing collision energy, y/s, the resulting fireballs tend 
to freeze out at decreasing values of the chemical potential. This is easy to understand, since the 

9 In the chiral limit, the smallest value of )1r is zero, and is achieved exactly at T c for /i# = 0. 
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amount of generated entropy (heat) grows with y/s while the net baryon number is limited by that 
number in the initial nuclei. 

The information about the location of the freezeout point for given experimental conditions 
is obtained by measuring the ratios of particle yields (e.g., baryons or antibaryons to pions), and 
fitting to a statistical model with T and }Xb as parameters. Such fits are amazingly good [fiQ], and 
the resulting points for different experiments are shown in Fig. Q 

As with any critical point, measurement of fluctuations can be used to determine when the 
system is in the vicinity of the critical point. By measuring variables sensitive to the proximity 
of the critical point as a function of monotonically increasing y/s of the collision, and observing 
non-monotonic dependence, one discovers the critical point [p9|]. The values of T\Ib corresponding 
to the freezeout at such a value of yfs give the coordinates of the critical point. 

As a concluding remark, it should be pointed out that the physics of the critical point is uni- 
versal (as far as slow and long distance phenomena are concerned), which allows to define certain 
experimental signatures independently of microscopic description. However, the position of the 
critical point on the phase diagram is determined by microscopic physics and is not universal at all. 
This obviously makes it very difficult to predict the coordinates of the critical point reliably as it 
is evident in the scatter of predictions in Fig. f|. On the other hand, the same fact should turn the 
knowledge of the position of the critical point obtained on the lattice, or in the experiment, into a 
powerful constraint on possible models of QCD thermodynamics. 
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